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Abstract 

We analyze the specific heat variations as a function of an external magnetic field 
of a simple model of superlattice that includes (i) in-plane ferromagnetic exchange, 
(ii) interplane ferromagnetic exchange, (iii) dipolar interactions, (iv) magnetocristalline 
anisotropy. The calculations are carried out at the spin wave level. The interplay be- 
tween the existence of a canting transition and the anisotropy of the structure generate 
non trivial behavior for the spin wave contribution to the low temperature specific heat as 
a function of an external magnetic field when dipolar interactions and magnetocristalline 
anisotropy are taken into account. 
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1 Introduction 



Magnetic multilayers made of alternating ferromagnetic films and non magnetic films (being the typ- 
ical thickness of these films several atomic monolayers) present interesting properties in the presence 
of an external magnetic field, as for instance the giant magnetoresistance effect found in [jij and also 
in other multilayers [0]. Up to now, many different measurements have been performed Hj|], but 
to our knowledge, the specific heat of magnetic multilayers has never been studied due to the ex- 
perimental difficulty of measuring a mass as little as a \ig. Recent developments in nanocalorimetry 
allow now to achieve the necessary resolution for such a measurement ||, ||. The aim of this arti- 
cle is to present a theoretical investigation of the low temperature specific heat of bulk multilayers. 
Two contributions to the specific heat are present: first, the contribution of paramagnetic layers and 
second, a small contribution of spin waves. We analyze here the variations of the low temperature 
spin wave contribution as a function of the strength of magnetic field applied perpendicular to the 
planes or parallel to the planes. This contribution could possibly be measured via a high accuracy 
specific heat measurement. Our goal is not to obtain quantitative accurate results, but rather to 
give a scenario for the behavior of the spin wave excitations and the consequences for the spin wave 
contribution to the specific heat. 

We thus work with a simple model including (i) in-plane ferromagnetic exchange of the order of 
J|| ~ 100 K. (ii) interplane antiferromagnetic exchange of the order of J± ~ 1 K. (iii) in-plane dipolar 
interactions with a strength of the order ~ 1 K. (iv) magnetocristalline anisotropy of the order 
K ~ 0.5 K. The magnetic moments are taken to be localized. This is only an approximation, but 
widely used in the literature (see for instance ||] ) . For simplicity, we study the magnetic multilayers 
as a bulk system, and we use periodic boundary conditions in each direction. In particular, we do 
not take into account the surface spin wave modes (see reference |8| and references therein). Even at 
this level of approximation, we find a non trivial behavior for the spin wave part of the specific heat, 
associated to the transition from a canted spin configuration to an aligned spin configuration as the 
strength of the magnetic field increases. This transition was already noticed (see || and references 
therein). We also find that, when dipolar interactions and anisotropy are taken into account, the 
anisotropy of the structure implies an anisotropy of the specific heat with respect to the magnetic field 
direction. Namely, the specific heat with a perpendicular magnetic field is larger than the specific 
heat with an in-plane magnetic field. 

2 Description of the model 

Let us now describe the model and fix the notations. The x and y axis are chosen parallel to the 
layers and the z axis is perpendicular to the layers. Since we are interested in bulk properties, we 
use periodic boundary conditions in each direction. Moreover, the spins lye on a cubic lattice. The 
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Hamiltonian reads 
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where the exchange tensor A is 
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the dipolar tensor being 



d a)/ a(x,x') 



1 



3 (8a,/3 ~ 3u Q ,(x,x').U/3(x,X / )) . 



(5) 



x' — X 



e x ,e y ,e z are the unit vectors in the x, y and z directions, and u(x, x') is the unit vector in the direction 



one should remember that h is related to the "true" magnetic field H by h = MH, with M the net 
magnetic moment. As far as dipolar interactions are concerned, we take into account only in-plane 
dipolar interactions. This is justified by the structure of magnetic multilayers: the ferromagnetic 
layers are separated by paramagnetic layers, so that it is reasonable, in a first approximation, to 
assume that dipolar interactions are long ranged inside a layer and that two spins in a different layer 
do not experience dipolar interactions. 

The term (|j) is a magnetocristalline anisotropy term, modeled via a short range pseudo dipolar 
interaction, as proposed in 0. We choose the anisotropy to be small enough so that in a zero 
magnetic field the magnetization is in the planes. This amounts to consider only the case of magnetic 
multilayers for which the thickness of the ferromagnetic layers is large enough [fTcl , typically larger 



3 Spin wave spectrum in a zero magnetic field without dipolar 
interactions and anisotropy 

This section is devoted to the zero magnetic field case. We determine the spin wave spectrum in 
a zero magnetic field and show that, due to the three dimensionality of the system, no divergences 
appear at the spin wave level (unlike the one dimensional Heisenberg chain case). 

3.1 Spin wave spectrum 

We first analyze the zero magnetic field spin wave spectrum and the specific heat of the "minimal" 
model including only the ferromagnetic Jn and antiferromagnetic Jj_ exchange constants. This model 
in a zero magnetic field can be solved at the spin wave level using the Holstein-Primakov bosons 
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is an external magnetic field. In order to make contact with the literature, 
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method [11]. Introducing two types of bosons a and b for each type of layer, keeping as usually only 
the 1/s leading quantum corrections in the semi classical (large s) approach and Fourier transforming 
the Hamiltonian leads to 



with 



^4q = J\\ (2 — cos q x — cos q y ) + J± 
Bq = J±cosq z . 

The Hamiltonian is readily diagonalized via the following Bogoliubov transformation 

c q = cosh y? q a q — sinh V2 q &l q 

dq = COsh <^ q 6 q — Smh V3qtlq, 



with tanh (2ip c 



-Bq/Aq, leading to the following spin-wave spectrum lu„ = A^ — B^ 



Jii (2 — cos q x — cos q y ) + J± J — J± cos 
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If J|| = 0, we recover the spin-wave spectrum w q = J± \ siaq z \ of the Heisenberg chain |1^] 13] and 
if Jj_ = 0, we recover the w q = Jn(2 — cosq x — cosq y ) spectrum of the two-dimensional ferromagnet 

& 

3.2 Stability of the classical alternate ground state 

Before going further, we would like to show that the antiferromagnetic perpendicular exchange J± 
does not destabilize the alternate ground state, unlike the one-dimensional case where the Neel state 
is destabilized by the first quantum corrections. The stability in the multilayer case is due to the three 
dimensionality of the structure. We calculate the first quantum corrections to the magnetization. To 
do so, we notice that 
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(12) 



leading to the average quantum correction AM to the magnetization per spin in a layer: 
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where BZ denotes an integral over the first Brillouin zone. This integral is infrared convergent: 
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where we have expanded the and l? q in the small q limit. A dimensional analysis shows that 
AM ~ J±/Jn. In the limit Jn — > 0, we recover the spin-wave divergence of the Heisenberg chain [12|. 
In the multilayers, Jy is two orders of magnitude larger than J±: our spin wave approach is thus 
consistent from the point of view of the absence of divergences of the spin wave theory. On the other 
hand, the three-dimensionality of the structure also avoids problems related to the Mermin- Wagner 
theorem in two dimensions at a finite temperature (absence of spontaneous order in two dimension 
for continuous spins). We have shown the consistency of our spin wave approach only in the case 
of a zero magnetic field and without dipolar interactions. However, the spin wave approach is still 
consistent in the presence of a magnetic field and with dipolar interactions. 



4 Spin wave theory in a perpendicular magnetic field 

4.1 Classical ground state without dipolar interactions and anisotropy 

We are first going to examine the spin wave theory when an external magnetic field is switched on, but 
in the absence of dipolar interaction and anisotropy. The effect of dipolar and anisotropy interactions 
will be discussed later on. At the classical level, the zero temperature ground state corresponds to the 
condition that each spin is exactly aligned in its local field. Because of the translational invariance, 
each spin in a given layer points in the same direction (at the classical level and at zero temperature). 
We denote by the spin vector in layers of type (a), and the spin vector in (b) type layers: 
the spins can be canted in a layer by layer fashion. Layers (a) and (b) are alternating along the z 
direction. The local field h^'^ in layers (a) and (b) respectively is given by 

hW = h + 2J||n( a ) - J ± n( 6 ) (15) 
= h + 2J||nW- JjnW, (16) 

where h is the external magnetic field. We denote by u the unit vector in the direction of the external 
magnetic field: h = \h\u. We are looking for a classical ground state such as n^ a \ and u are 
coplanar. As we will see below, such a classical spin configuration exists. We denote by a the angle 
between the vectors u and n( a ): n( a ) = ucosa + vsina, where v is a vector orthogonal to u. The 
angle a between the magnetic field and the spin direction is called the canting angle in what follows. 
The local field h^-* experienced by a spin in a type (a) layer is easily expressed as 

h (a) = u (h + (2 Jy - Jl) cos a) + v (2J\\ + J_l) sin a = |h (a) | (u cos a + v sin a) . (17) 

The last equality simply means that, at the classical level, the spins are aligned with their local field. 
We can thus determine the direction of the n( a ) vector: 

cos a = h/2J± iih<2J± . » 

a = if h > 2J±. ' 

At the classical level, the spins are thus exactly aligned with the magnetic field provided the intensity 
of the magnetic field is larger than 2J± and in the absence of dipolar interactions and anisotropy. 
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4.2 Spin waves spectrum without dipolar interactions and anisotropy 

In the absence of dipolar interactions and anisotropy, the classical canted configuration of spins is 
given by (|l8|). The in-plane spin wave degrees of freedom decouple from the perpendicular degrees 
of freedom. The contribution of the in-plane spin wave degrees of freedom is a trivial ferromagnetic 
term, which will be added at the end of the calculation. For the moment we consider only spin 
wave excitations along the z direction, that is, we are going to solve the Heisenberg chain with an 
antiferromagnetic exchange J± in a magnetic field at the spin wave level. In order to obtain the 
appropriate excitations, we first need to express the spin operators a ±)Z in terms of the rotated spin 
operators a R ' z corresponding to a spin quantization along the direction of the classical ground state. 
The rotated spin states are 

\+)r = cos(|)|+)-*sin(|)|-> (19) 
\-) R = cos(!)|->-;sin(!)|+>, (20) 

where the rotated spin states |±)# are deduced from the |±) states by a rotation of angle a around 
the x axis. This allows us to express the spin operators in the laboratory framework in terms of the 
spin operators in the local framework: 

<t + = -^smaa z R + cos 2 (a/2)a^+sm 2 (a/2)a R (21) 

a~ = ^smaa R + cos 2 (a/2)a R + sin 2 (a/2) (22) 
a z = cos aa z R — i sin air| + i sin aa R . (23) 

The expressions should be plugged into the Hamiltonian of the Heisenberg chain in a magnetic field. 
The spin wave approximation then consists in using boson operators related to the local spin operators 
as follows: a R = a, a R = a + and a z R = 1 — 2a + a. Technically, the effect of the rotation to the local 
framework is to eliminate linear terms in boson operators. We are thus left with a hamiltonian that 
is quadratic in terms of Holstein-Primakov bosons. Considering only two consecutive spins with 
boson operators a and 6, and after straightforward calculations, we obtain the following spin wave 
Hamiltonian 

H (SW ) = J± ( a+a + b+h ^ + J± gin 2 a ( a+6 + + ah ) + l i L cos 2 a ( a6 + + a + fe ) . ( 2 4) 

After Fourier transforming the chain Hamiltonian and performing a Bogoliubov rotation, we deduce 
the following spin wave spectrum uj 2 z = A 2 z — B 2 z with 

A^> = Jj_ (l + cos 2 a cos (25) 
B^~' = J_i_ sin 2 a cos q z . (26) 
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This expression is defined on the unfolded Brillouin zone [— tt, tt]. However, one should remember 
that, because of the doubling of the unit cell, the "true" Brillouin zone is [— vr/2, 7r/2]. In order 
to come back from the unfolded Brillouin zone [— vr, vr] to the "true" Brillouin zone [— vr/2, vr/2], one 
should consider two spin wave modes for a given wave vector. Notice also that the unfolded dispersion 
relation is not unique: starting from two distinct spin wave modes in the "true" Brillouin zone, there 
are two ways to unfold the dispersion relation to the [— vr, vr] interval. The spectrum defined by 
(25) and (26) corresponds only to one of the unfolding procedures. The other unfolding procedure 



corresponds to changing q z into q z + vr in the unfolded Brillouin zone. 

We notice that, in a zero magnetic field, there are two Goldstone modes. As the magnetic field 
is switched on, one of these Goldstone modes acquires a gap while the other one remains massless 
(rotations around the magnetic field). The dispersion relation defined by ( p5|) and (|2^ ) is thus 
consistent with the fact that, at the classical level in the presence of a non zero magnetic field the 
classical ground state is defined up to a rotation around the magnetic field, thus leading a single 
gapless mode at the spin wave level. 



We have rederived the spin wave spectrum defined by (25) and (26) using a method based on 
the semiclassical equations of motion, which lead exactly to the same result. We do not present here 
these calculations since they do not lead to new physical insights. 

In order to come back to the multilayer problem, we have to take into account the existence of 
in-plane ferromagnetism, leading to an additional A term that should be added to (p5|): 



^11) = J|| (2 - cos q x - cos q y ) . (27) 

Now, if h > 2Jj_, the spins are aligned with the magnetic field. The dispersion relation is 
straigthforwardly found to be 

iO q = h — J± + J± cos q z + J|| (2 — cos q x — cos q y ) . (28) 

This leads to a gap A = h — 2J±. 

We have thus obtained an analytic expression for the spin wave spectrum in the absence of dipolar 
interactions and anisotropy. The evolution of the dispersion relation for q x = q y = as a function of 
q z and the magnetic field is plotted on figure |l|. 

We emphasize that this result does not depend on the direction of the magnetic field: since the 
Hamiltonian without magnetic field is spin rotational invariant, the spin wave spectrum should be 
the same whatever the direction of the magnetic field. 

4.3 Specific heat without dipolar interactions and anisotropy 

The specific heat per spin is obtained from the dispersion relation via 

' (29) 



4T 2 J M s inh 2 (/3w q /2) ' 
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where j3 = 1/T is the inverse temperature. We calculated numerically the spin wave specific heat as 
a function of a magnetic field for a low temperature T = 3K. The result is plotted on the curve (a) 
of figure 0. We observe the existence of two regimes: (i) h < h c = 2J±: the specific heat is slowly 
decaying as the magnetic field is increased. We attribute this slow decay to the fact that the system 
remains gapless in the presence of the magnetic field. If h = 0, there are two Goldstone modes. One 
of these Goldstone modes becomes massive as the magnetic increases, leading to a decrease of the 
low energy density of states and thus a slow decrease in the specific heat, (ii) h > h c : the system 
acquires a gap A = h — 2J±, leading to a decay to zero of the specific heat as the magnetic field 
increases, with a cross-over magnetic field h* = 2J± + T. 

4.4 Spin wave spectrum in the presence of in— plane dipolar interactions 

We now assume the existence of in-plane dipolar interactions. This assumption is motivated by the 
fact that the ferromagnetic layers are separated by paramagnetic layers, and it is thus reasonable to 
suppose that only spins in the same layer are coupled by dipolar interactions. 

The first step is to calculate the zero temperature classical configuration of spins. This can be 
easily done by minimizing the classical energy with respect to the canting angle and we find 

h 

cosq = (30) 

2 J ± + 3vr J d y ' 

in the regime where the spins are canted. Physically, we expect that, at the classical level, the 

effect of dipolar interactions is to open the canted angle since dipolar interactions play in favor of 



attracting the spins in the planes. This is indeed the case, as shown by equation (30). Notice that in 
the evaluation of the dipolar interaction energy, we have replaced the discrete summation over the 
lattice by an integral. We will use this approximation throughout this paper to treat the effects of 
long range in-plane dipolar interactions. 

Next, we need to calculate the spin wave spectrum in the presence of dipolar interactions. The 
calculations are lengthly but straightforward so that we just give here the results. The idea consists in 
eliminating linear Holstein-Primakov bosons by a rotation to the local framework. Keeping only the 
quadratic terms in boson operators leads to the following results in the canted phase. The additional 
contribution to the A q term due to dipolar interactions is 

Af*) = ^ cos 2 a + nJ d (l (q l{ ) - l) + ^^o(<?||) (Ssin 2 a - 2) + ^ sin 2 acos (2^)/ 2 ( 9 ||) > 

(31) 

where 4> q is the angle between the vector (q x ,Qy) and the x axis, and I n is expressed in terms of 
Bessel functions as 

7 n(?) = J i ^Uqr). (32) 
The dipolar contribution to the B q terms is 

B (di P ) = _3TrJd ( sin 2 aJo ( g|| )_ cos ( 2 g )i 2 ( g|| ) ((1 + cos 2 a) -2i cos a)) . (33) 
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As expected, if q x = q y = 0, |Aq^| = \Bq | since the q x = q y = mode in a single ferromagnetic 



layer with dipolar interactions is gapless. Notice that the contribution (33) is in general a complex 
number. However, this is not a problem: using a phase transformation on the Holstein-Primakov 
bosons, we show that the spin wave spectrum is given by uj 2 = (A 1 " 1 ) 2 — \B^ ot \ 2 , where the subscript 
"tot" refers to the total A and B terms (the contributions of ferromagnetic in-plane exchange terms 
+ antiferromagnetic coupling terms + dipolar interaction terms have been added). Specializing the 
dispersion relation to the q x = q y = case, and in the absence of a magnetic field (a = 7r/2), 
we observe the existence of two Goldstone modes in the absence of dipolar interactions but in the 
presence of dipolar interactions, only one of these modes remains massless whereas the other one 
acquires a gap. This remaining massless mode corresponds to rotations of the classical spins inside 
the planes. In the presence of a magnetic field, this massless mode remains since it also corresponds 
to the rotation of the spins around the magnetic field (perpendicular to the planes). 

Notice also that the terms proportional to ^(^ll) m (pl|) and (|33| ) are expected to lead to a very 
small (and neglectable) contribution to the specific heat: the main contribution to the specific heat 
comes from the lowest energy states, namely the states next to q x = q y = 0. Along this line, the 
terms proportional to li vanish since ^(0) = 0. In the rest of the Brillouin zone the contribution of 
the terms involving I2 remains very small compared to the in-plane ferromagnetic exchange term, 
since Jy is two orders of magnitude larger than J^. In the present case of a perpendicular magnetic 
field, the terms proportional to I2 can be treated exactly since they are the same for all types of 
layers (invariance under the a — > —a transformation). As we shall see below, this is not the case 
for a magnetic field parallel to the layers. However, we will argue that one can safely neglect these 
terms since they lead to a neglectable contribution to the specific heat. 

The dispersion relation on the line q x = q y = as a function of q z and the magnetic field is 
plotted on figure ||. 

4.5 Magneto cristalline anisotropy as a pseudo dipolar interaction 

We now take into account an anisotropy term (Q), under the form of a pseudo dipolar in-plane 
interaction. We just give the results. One should first calculate the canting angle in the presence of 
such an interaction. After minimizing the classical energy, we find easily 

cos a = ; . (34) 

2J ± + 3vr J d /2 -3K K J 

This equation is consistent with the fact that the pseudo-dipolar interactions play in favor of aligning 
the spins perpendicular to the layers, and thus reduce the canting angle. As far as the spin wave 
spectrum is concerned, the first effect of these pseudo-dipolar interactions is to renormalize Jm: 
J\\ ~~ * J\\ ~ Since K is much smaller than Jm, this effect is very small. The contribution of the 
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anisotropy to A is found to be 



A q K ^ = —3K H — — ^cos q x + (cos 2 a) cos q y J , (35) 
and the contribution to B is 

B^ = — (cos q x - (cos 2 q) cos q y j . (36) 

As in the case of dipolar interactions, I-Aq^I = l-^o I ^ Qx = Qy = since the q x = q y = mode 
is a Goldstone mode of a ferromagnetic layer with pseudo-dipolar interactions. Figure |] shows the 
dispersion relation along the line q x = q y = as a function of q z and the magnetic field. Qualitatively, 
adding an anisotropy term does not change the main features of the dispersion relation (namely the 
existence of a single massless Goldstone mode). 

4.6 Specific heat in the presence of in— plane dipolar interactions and magneto 
cristalline anisotropy 

We are now in position to calculate the variations of the low temperature specific heat as a function 
of the intensity of a magnetic field applied perpendicular to the layers, within the spin wave approx- 
imation for the Hamiltonian (§)@@- These variations are plotted on the curve (b) of figure |. We 
observe two regimes: (i) for h < h c = 2J± + 37rJ^ — 3K, we observe that the specific heat increases as 
the magnetic field increases. We attribute this behavior to the fact that increasing the magnetic field 
induces an increase of the low energy density of states. This can be seen as follows. If q x = q y = 0, 
the dispersion relation around the Goldstone mode q z = tt — Aq z (in the unfolded Brillouin zone) is 
of the form 

uj qx=0 ,q y=0t q z=n -Aq z =Xs'maAq z . (37) 
with A independent on the magnetic field: 

^(^-¥)- <*> 

The density of states in the vicinity of the massless mode is thus an increasing function of the 
magnetic field, leading to an increase of the specific heat as the magnetic field increases. 

(ii) h > h c : the gap is A = h — h c thus leading to a decrease to zero of the specific heat as the 
magnetic field increases with a cross-over magnetic field h* = T + h c = T + 2Jj_ + 37rJ^/2 — 3K. 

Finally, we can compare this result with the specific heat in the absence of dipolar interactions. 
In a zero magnetic field, we expect the specific heat without dipolar interactions and anisotropy 
to be larger than the specific heat with dipolar interactions and anisotropy since there is a single 
massless mode in the former case and two massless modes in the latter. This is indeed the case, as 
shown on figure ||. Now, consider the aligned state (h > h c ). The gap is reduced when dipolar and 
pseudodipolar interactions are switched on, so that switching on dipolar interactions and anisotropy 
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should increase the specific heat, which is what is observed on figure The existence of a crossing 
point between the specific heat curves with and without dipolar interactions and anisotropy on figure 
H is thus consistent with our qualitative analysis. 



5 Spin wave theory in a longitudinal magnetic field 

We now examine the case of a magnetic field applied parallel to the layers. The effect of the orientation 
of the magnetic field can only be seen when dipolar interactions and anisotropy are taken into account. 
Namely, in the absence of dipolar interactions and anisotropy, the spin wave spectrum is the one 



calculated in section 4.2. We are going to calculate the effects of a longitudinal magnetic field. To do 



so, we use the same technique as in section we first express the Hamiltonian in terms of the spin 
operators quantized along the local magnetization axis, which leads to a quadratic Hamiltonian at the 
spin wave level, that can be diagonalized via a Fourier transform and a Bogoliubov transformation. 
First, we treat the problem in the absence of dipolar interactions and anisotropy with the spins in 



the planes direction. Of course, we recover the results of section 4.2. In a second step, we discuss the 
effects of in-plane dipolar interactions and anisotropy. In order to avoid confusion with the previous 
calculations, we call ip the canting angle. 

5.1 Spin wave spectrum in the absence of dipolar interactions and anisotropy 

If w denotes the unit vector of coordinates w(sin ip, cos ip, 0), with (p the canting angle, then, the 
local spin states in a given layer are 

\+) R = J=(|+) +e -^|-)) (39) 
|->* 



v/2V 

^ (-e*|+> + |-)) , (40) 



where we have performed a rotation of tt/2 around the w axis. The case of the other type of layer 
is obtained by replacing ip by — (p in these equations. The expression of the spin operators a ±,z in 
terms of the local spin operators a R ' z is 



(7 



»r) ( 41 ) 
\ (f^R + *Z ~ ^r) ( 42 ) 



d» = - e ^a+-e-^* R . (43) 



Using the same technique as in section 4.2, we express the Hamiltonian of a single a-b link in terms 
of the a and b bosons (one type of boson for each direction of the classical ground state). In order to 
get rid of phase factors, one has to define new boson operators 5 = aexp (iip) and b = ftexp (— itp). 
We then recover the spin wave Hamiltonian (24) in terms of the d and b bosons. 
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5.2 Spin wave contribution of in— plane dipolar interactions 

The first step is to determine the classical spin configurations. In the presence of dipolar interactions, 
the classical spins are aligned in the planes, with a canting angle 

cos</? = ^-. (44) 

Notice that the dipolar interaction coupling does not come into account in the expression of the 
canting angle. This is due to the fact that the spins are in the planes, and therefore the classical 
energy of dipolar interactions is a constant, independent on the canting angle <p. This is already a 
first difference between the cases of a perpendicular and longitudinal magnetic fields. 

The second step is to include in-plane dipolar interactions at the spin wave level. We do not give 
the details of the calculations. The contribution of dipolar interactions to A for a given layer is 

A<&> = ttJ, (j ( g ,|) - l) + ^Jo(g||) - ^ cos (2(0, - ^))I 2 (<?||), (45) 
and the contribution to B is 

Bf p) = ^ (loin) + cos W 9 + V)H 2 (<7||)) • (46) 

These calculations have been carried out in a given layer. The main problem regarding these ex- 
pressions is that they are not invariant under the (f — > —ip transformation. In other words, the 
dipolar terms induce a hybridization between the two Goldstone modes in the [— 7r/2,7r/2] Brillouin 
zone. However the strength of the hybridization is proportional to l2(q\\). As mentioned earlier, we 
expect the terms proportional to Ii to lead to a really small contribution to the specific heat since 
i) I2 vanishes along the line q x = q y = 0, where the lowest energy modes are located ii) outside 
the q x = q y = line, these terms are dominated by the ferromagnetic exchange interaction which 
is two orders of magnitude larger than the dipolar coupling. We can thus safely neglect the terms 
proportional to I2 in our specific heat calculation. 

Finally, in the presence of dipolar interactions, the results in a perpendicular and longitudinal 
magnetic field should be the same in the limit of vanishing fields. This is indeed the case. In order 
to get the same result for a = ip = vr/2, one should make a ir shift on the z component of the wave 
vector. This ir shift is allowed since, as mentioned earlier, there are two possibilities of unfolding the 
"true" Brillouin zone [— 7r/2,7r/2] to the extended Brillouin zone [— ir,Tv]. These two possibilities are 
precisely related to each other by a tt shift on the z component of the wave vector. 

Before discussing the results (03) and (|46|), we would like to include the effects of an anisotropy. 



5.3 Spin wave contribution of the magneto cristalline anisotropy as a pseudo 
dipolar interaction 

The magnitude of the anisotropy is chosen to be small enough so that the classical configuration 
of spins is aligned in the plane, whatever the value of the longitudinal magnetic field. Under this 
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assumption, the canting angle is given by (fl||) and is independent on the strength K of the anisotropy. 
The anisotropy term can straightforwardly be incorporated in the spin wave approach. The results 
are as follows. The first effect of the anisotropy is a renormalization of Jm: Jii — ► Jii — K. This effect 
is very small since J11 is two orders of magnitude larger than if. The contribution of the anisotropy 
to the A term is 

= -3K + — (sin 2 tp cos q x + cos 2 </? cos q^j . (47) 
The contribution of the anisotropy to the B term is 

B^ = — — (sin 2 ip cos q x + cos 2 (p cos q z J . (48) 

The dispersion relation in the presence of an anisotropy and dipolar interactions is plotted on figure 
[| in the extended Brilouin zone for q x = q y = 0. In a zero magnetic field, we observe the presence 
of a massless mode associated to the rotation of the spins inside the planes. As soon as a magnetic 
field is switched on, this mode acquires a gap, due to the fact that the magnetic field breaks the 
rotational invariance of the ground state with respect to in-plane rotations. Another effect is that 
if h = 2J±, we observe the existence of a massless mode. The existence of this gapless mode can 
be understood at the classical level as follows: consider a magnetic field strictly equal to 2J±. It is 
then possible to calculate the energy of a classical spin wave of amplitude 9 around the aligned state, 
and expand the energy in powers of the amplitude 9 of the classical spin wave. One then get the 
result that the leading term in this expansion is of order 9 4 . More precisely, the energy per spin of 
an antiferromagnetic classical spin wave of amplitude 9 around the aligned state (h > 2J±) is 

E(9) = -^ cosd + ±L cos (29) = -\{h- 2J1)(1 - 9 2 ) - - 8J ± )0 4 + ... (49) 

If h = 2J± the leading term of this expansion is of order 9 4 , thus leading to a flat mode at the 
harmonic order. It is thus not surprising that at the quantum spin wave level one finds a gapless 
mode right at the canting transition. 

5.4 Specific heat in a longitudinal magnetic field 

The specific heat in a longitudinal magnetic field is plotted on the curve (c) of figure We observe 
the same type of variations as in the perpendicular magnetic field case, except that the specific heat 
with an in-plane magnetic field is smaller than the specific heat without dipolar interactions and 
anisotropy. This can be explained qualitatively as follows. In the aligned state (h > 2J±), the gap is 
given by 

A 2 = (h - 2J ± ) 2 + (h - 2J ± ) - 3^, (50) 

where the dispersion relation has been specialized to the q x = q y = 0, q z = it case. As a consequence, 
the gap with an in-plane magnetic field and dipolar interactions and anisotropy is larger than in the 
absence of dipolar interactions and anisotropy. We thus expect from this argument a lower specific 
heat, which is indeed the case as observed on figure |2[ 
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6 Conclusion and discussion 



Let us now summarize our results. We have calculated the spin wave contribution to the specific heat 
in three different cases: (a) without dipolar interactions and anisotropy. At this level of approxima- 
tion, the specific heat does not depend on the direction of the magnetic field since the Hamiltonian 
without a magnetic field is spin rotational invariant, (b) with dipolar interactions and anisotropy for 
a perpendicular applied magnetic field, (c) with dipolar interactions and anisotropy for a longitudinal 
applied magnetic field. In case (b), the location of the canting transition depends explicitly on the 
strength of dipolar interactions and anisotropy. In case (c), the location of the canting transition is 
independent on the strength of dipolar interactions and anisotropy. We observe a maximum of the 
specific heat at the canting transition in cases (b) and (c). Another conclusion is that, due to the 
anisotropy of the system, we predict an anisotropy in the specific heat. Namely, the specific heat in 
a perpendicular magnetic field is larger than the specific heat in an in-plane magnetic field. 

However, there is also a contribution to the specific heat due to non magnetic interlayers. As 
only a few percentage of their conduction electrons is polarized by the exchange coupling |l(| the 
interlayers should behave nearly as a bulk paramagnet. This paramagnetic contribution per spin has 
then the form 

Jpara) _ (ffi) (t --,\ 

v "cosh 2 {fihy 1 j 

where f3 is the inverse temperature. This paramagnetic contribution is much larger than the spin wave 
contribution which is a source of experimental difficulty. Clearly, the sample should be processed in 
such a way as to minimize the number of paramagnetic atoms. 
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Figure captions 



Figure 1: 

Dispersion relation as a function of the magnetic field. The dispersion relation is plotted along the 
line q x = q y = as a function of the strength of the magnetic field in the extended Brillouin zone 
[— vr,7r]. There is no dipolar interaction nor exchange interaction. We took J± = 1. 



Figure 2: 

Specific heat per spin in the different cases considered in this article, with J± = IK, Jy = WOK, and 
as a function of the external magnetic field. The curve (a) corresponds to the specific heat in the 
absence of dipolar interactions and anisotropy. (b) corresponds to a perpendicular applied magnetic 
field, with Jd = IK and K = 0.5K. (c) corresponds to a longitudinal magnetic field with = IK 
and J ± = 0.5K. 



Figure 3: 

Dispersion relation as a function of the intensity of a perpendicular magnetic field. The dispersion 
relation is plotted along the line q x = q y = as a function of the strength of the magnetic field in 
the extended Brillouin zone [— w, 7r]. We took Jj_ = 1, = 1 and K = 0. Notice the differences with 
figure ffl. 

Figure 4: 

Dispersion relation as a function of the intensity of a perpendicular magnetic field. The dispersion 
relation is plotted along the line q x = q y = as a function of the strength of the magnetic field in 
the extended Brillouin zone [— 7r,7r]. We took J± = 1, = 1 and K = 0.5. The variations of the 
dispersion relation are very similar to the ones of figure ||| where K was set to zero. 

Figure 5: 

Dispersion relation as a function of the intensity of a longitudinal magnetic field. The dispersion 
relation is plotted along the line q x = q y = as a function of the strength of the magnetic field in 
the extended Brillouin zone [— vr,7r]. We took J± = 1, = 1 and K = 0.5. Notice the differences 
with figure |]. 
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